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1. Introduction 

Relatively few hadronic states are stable to strong decays (i.e. via QCD with degenerate u and 
d quarks). Among the mesons, we have [p: 

Stable n K r\ D D s B B s B c D* B* B* s D, (0+ ) B, (0+ ) 

T < 1 MeV T\' D* \jf(lS) y/(2S) X\ Xi T(15) T(25) T(35) 

r<10MeV co (j> XoX(3S72) 

r > 10 MeV p /o ao &i &i ai /2 /i ^2) etc., inc r\ c . 

The mass of an unstable state is usually defined as the energy corresponding to a 90° phase 
shift. This definition seems to accord with simple mass formulae: For example 

• p(776) and to (783) are close in mass despite having widths of 150 and 8 MeV respectively. 

• The baryon decuplet (A(1232), 1(1385), 2(1530), £2(1672)) is roughly equally spaced in 
mass despite having widths of (120, 37, 9, 0) MeV respectively. 

(Note that defining the mass as the real part of the pole will cause a downward shift of masses for 
wider states, eg. 22 MeV less [g| for the A(1232) pole, and this prescription will fit the equal mass 
rule less well.) 

So, on the one hand, unstable particles seem to fit in well with stable ones; on the other hand, 
the presence of open decay channels will have an influence in lattice studies. 

Some of the motivations to study hadronic decays on the lattice are: 

• To determine properties of exotic states (glueball; hybrid meson; multi-quark) to guide ex- 
periment. 

• To understand the nature of states: whether meson-meson or quark-antiquark in structure. 

• To understand the impact of decay on the mass of a state. 

• Hadronic decays are now accessible in full QCD: for example po — ^ TTi 7T_ i at rest requires 

m(n) 0.5 

< 



m(p) Ai I 4k 2 ~ 

so for m{%)L = 5 this implies m(%)jm[p) < 0.32, but note that pi — ► %\Xq , where suffix 
labels momentum in units of 2n/L, only needs m(7t)/m(p) = 0.44, 

2. Decays in Euclidean Time 

NO GO. At large spatial volume, the two-body continuum masks any resonance state. The 
extraction of the spectral function from the correlator C(t) is ill-posed unless a model is made Q], 
since the low energy continuum dominates at large t. 
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GO. For finite spatial volume (L 3 ), the two-body continuum is discrete and Liischer showed [Q, 
|6[ 0] how to use the small energy shifts with L of these two-body levels to extract the elastic 
scattering phase shifts. The phase shifts then determine the resonance mass and width, see ref. [||] 
for a review. Thus a relatively broad resonance such as the p appears as a distortion of the % n Tl- n 
energy levels where pion momentum q = 2%n/L. 

The effect can be visualised as arising from the relatively larger amplitude for interaction 
between two hadrons in increasingly smaller spatial volumes: 



o 



o 




3. Lattice evaluation 

Liischer's method is applicable in principle to study phase shifts below inelastic thresholds. 
But problems arise in practical implementation: 

• Accurate measurement of a small energy difference (of order L~ 3 in general, though bigger 
at resonance) is needed. 

• Measurement of a matrix of correlations between two-body and one-body operators will 
be needed to get accurate energy values and these correlators are less straightforward to 
evaluate. 

• Measurements are needed at several spatial sizes (L) with dynamical fermions, with L big 
enough to hold two hadrons. 

• Higher energy levels in a given channel are particularly hard to determine (eg. %\ 7T_ i com- 
pared to KqKq) since excited state energies are always difficult to extract with precision from 
temporal correlators, since the ground state contribution dominates at larger t. 

Are there any short cuts? 

4. Lattice evaluation: tricks 



4.1 7T+7T+ 

One of the first cases that has been explored is the n + 7l + interaction at low energy. This is 
favourable since there are fewer quark diagrams to evaluate as there is no annihilation diagram, 
i.e no qq channel and, because of this, quenched evaluation is a useful approximation. This has 
been much studied in the last decades. Indeed recently there have also been dynamical quark 
evaluations [^, 10] of the n + 7t + scattering near threshold. 

The trick that has been evolved to improve precision in these studies is that of using an operator 
that creates two pions at spatial separation x. Then the nature of the expected spatial wave function 



008/3 



Hadronic Decays 



C. Michael 



versus x allows to determine the scattering length (phase shift 5/=2 near threshold) more accurately 
than from the energy determined by the £-correlation alone [pQ. 

This n + 7i + case does not involve decay, of course, just an evaluation of the hadronic interac- 
tion strength. To explore decays, one must study the transition between two meson and one meson 
operators. 



4.2 Hadronic transitions 

Consider a lattice study of the off-diagonal correlator: from a p meson to xx. Diagrammati- 
cally: 

p — > XX 

X 

t T 

Now to evaluate this contribution, since the intermediate point marked X at time t is not ob- 
served on a lattice, it must be summed over. 

T 

Y^e^ m[p), xe^ m{ - n7l)(T ^ , ' ) -^Ae~ m(p)T -Be~ m{n7l)T 

The problem that arises is that excited states of either the p or of xx will contribute a similar 
behaviour: 



-m'(p)t xe -m(nn)(T-t) 



t=0 



. Ce -m'{p)T _ De -m{nn)T 



J^ e -" 1 ^ t X e- in ' ( - 7t7C ^ T -^ 



Hence there is no way in principle to remove any excited state contamination unless m(p) « 
m(xx) when the ground state piece sums to xTe~ mT while excited states only behave as e~ mT . 

Thus for on-shell transitions on the lattice, it is possible to extract the hadronic transition 
amplitude directly [12, 13, 14]. The key signal is to observe a linear dependence of the lattice 
normalised transition amplitude on the temporal extent T. The slope of this transition then gives the 
lattice amplitude x that can be related to the transition amplitude with conventional normalisation 
of states. 



5. Lattice evaluation: p — >• %% 

This transition can be evaluated directly - for on shell transitions - by looking for the signal 
extensive in T. This has been explored [|l4|] for p decay to xx. In order to have an on-shell 
transition with the dynamical quark lattices then available, it was optimum to study the case of 
decay in flight: p\ — > X\Xq which is quite close to on-shell as illustrated in fig. 1. The diagrams 
illustrated in fig. 1 were evaluated using a stochastic time-plane source method. The generalisation 
of Luscher's method to decay in flight is given in JT^]. 
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id 




JtJI->JIJt 71 71 -> 71 71 



Figure 1: (a) Energy levels of p meson and ktz states. From UKQCD with Nf — 2 sea quarks with NP- 
improved fermions, m(n) / 'm(p) = 0.58, m(n)L = 4.6, a=0.1 1 fm. (b) Diagrams evaluated (stochastic time- 
plane source method) 




5 10 5 10 

t/a t/a 



Figure 2: Strength of transitions measured on the lattice. 



The signal obtained from the lattice study is illustrated in fig. 2 both for the normalised p — » 
7iK transition and for the normalised "box" contribution %% — > 7T7T which has a contribution from 
Tin —>■ p —> 71%. The approximate consistency of these two approaches is a useful cross-check. 

Assuming two close energy levels (with mixing amplitude x) mix to give an energy shift Am, 
then m(pi) moves down, while E{k\Kq) moves up. The lattice results are 

• Transition pi — ► K\Uq\ Signal xT . a Am = 0.022+y 7 

• Box 7Ti 7To — > 7Ti 7i:o: Signal (xJ) 2 . (consistent with above) 

• Shift of energy of pi spin parallel to momentum (which mixes with %\ 71q) to spin perpendic- 
ular to momentum (which does not). aAm = 0.026(7) 
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• Liischer shift (unbinding energy of K\Kq). aAin = 0.04(3). Note big error 

So the tricks give smaller errors than the direct determination of the Liischer mass shift. The 
tricks have some inherent systematic errors, however, and are most useful for a qualitative study. 

These mass shifts with the assumption of p dominance of the %% partial wave give a determi- 
nation of the coupling g pK7t . Thus one can determine this coupling constant from the lattice (where 
decay does not proceed) and compare with experiment: 



method 




m sea 


8 




Lattice xT 


s 


s 


1.401 


-47 
-23 


Lattice p shift 


s 


s 


1.561 


-21 
-13 




s 


11. d 


1.5 




K* — > Kk 


u, d/s 


u, d 


1.44 




p —>nn 


u, d 


u, d 


1.39 





The agreement is good, indicating that there is a relatively mild dependence of the coupling 
constant on the sea quark mass (which is higher in the lattice study than in experiment). This 
satisfactory confrontation between lattice and experiment, is encouraging for lattice exploration of 
cases where predictions need to be made since experimental data are not available. 

6. Hybrid meson decay 

One of the characteristic predictions of QCD is that there can be mesons in which the gluonic 
degrees of freedom are non-trivially excited. The simplest example is a hybrid meson with spin- 
exotic J PC = 1~ + which is a J PC combination not available to a qq state. 

6.1 Heavy quarks 




Figure 3: Strength of string de-excitation transitions measured on the lattice. 
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The cleanest environment in which to study such states on a lattice is in the limit of very heavy 
quarks - bb. This can be approximated by using static quarks and the gluonic excitation arises 
as an excited string state between these static quarks with non-trivial gluonic angular momentum. 
Lattice studies have long predicted the spectrum of such states. 

To guide experiment, however, it is important to know the expected decay mechanism and 
associated width. In the static quark limit, several symmetries can be used which imply Jl6| ] that the 
dominant decay will be string de-excitation (rather than string breaking). Lattice study [ |l6| ] shows 
that the dominant decay of Ht is string de-excitation to Xbfo- The transition that is considered is 
shown by the diagram on the left of fig. 3. The lattice data on the transition are illustrated in fig. 3: 
the transition is closest to on-shell for R « 0.2 fm. The width is predicted to be around 80 MeV. 

This estimate from first principles of the decay width is of significance in guiding experimental 
searches for such hybrid states. 

6.2 Light quarks 

For light quarks, evidence for contributions from two-body states (such as 7tb\) to the spin- 
exotic (J PC = 1~ + ) channel is seen in dynamical studies [|l7|] - this complicates the extraction of a 
spin-exotic hybrid meson. 

A recent exploratory study (quenched) of %a\ mixing with the J PC = 1~ + spin-exotic hybrid 
meson has given an estimate of the width for this decay [|TJ]. 



7. String breaking 

In quenched QCD, as the static quark and antiquark sources are pulled apart the potential 
energy rises as oR, where a is the string tension. In full QCD, a sea quark-antiquark can be 
created from the vacuum and it becomes energetically favourable to have a meson-antimeson pair 
beyond some value of R. This is the phenomenon of string breaking. It has long been realised that 



this can be studied as a mixing phenomenon [19] with channels QQ and Qq qQ. For static quarks 



at separation R, there will be a level crossing and associated mixing of V(R) and 2m(B). This 
mixing is the measure of string breaking [jfj|], [Q. This mixing energy is very hard to determine 
on a lattice: since the crossing occurs at relatively large R ps 1.25 fm where amplitudes are small. 
At this crossing point, for static quarks, there is a mixing transition which will be independent of 
lattice spatial size L for sufficiently large L since each state is localised (the heavy quarks are static). 
Employing "all the tricks in the book" a first estimate of this energy shift (mixing amplitude) 



has been obtained [ )21[ ] of 5 1(3) MeV. An illustration of the crossing region is in fig. 4. 

Using the adiabatic approximation, this energy gap can be used to evaluate amplitudes for 
excited T decay to BB. So it does indeed have some relevance to the topic of hadronic decays. 



8. Scalar Mesons 



8.1 Light quarks 

Since uu + dd, ss, glueball, and meson-meson components are all possible for flavour-singlet 
scalar mesons, this is a difficult area to study both on a lattice, and in interpreting experimental 
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Figure 4: String breaking (here a m 0.083 fm). 



data. For scalar mesons the lowest mass decay channels are nn (flavour singlet: /q) or r\n (flavour 
non-singlet: ao) and these decay channels are open in many dynamical lattice studies. The history 
of attempts to study the complex mixing between these different contributions is : 



++ Glueball decay — > Tin: quenched study [22, 23]. 



Glueball mixing with qq meson. This is a hadronic transition which is independent of L for 



large L, so a mixing energy can be quoted. It has been studied using quenched [ |24| | and 



dynamical lattices [12] 



A full lattice study is needed which includes glueball, qq and %% channels but the disconnected 
diagram for fo — > %% is very noisy in practice. To reduce the contribution from disconnected 
diagrams, one can study flavour non-singlet decays. The simplest case is ciQ—>rin and this has been 
explored in quenched studies which have an anomalous behaviour: since the v\ itself is unphysical 
(appearing as a double pole degenerate in mass with the pion). Rather than try to correct for this 
anomaly which gives a wrong sign to the ao correlator at larger t, it is preferable to use a ghost-free 
theory. With two flavours of sea quark (TV/ = 2), the ao — > f]7t transition is physical but it involves 
the evaluation of an additional disconnected diagram - as illustrated in fig. 5. This transition is then 
approximately on-shell and the transition strength can be evaluated, as shown in fig. 5. 

The results [ p5| ] are indeed as anticipated for the transition involving particles with no mo- 
mentum - showing an approximately linear rise. For the decay in flight the result is significantly 
different even though, for an S-wave transition, one would expect a transition amplitude indepen- 
dent of momentum. This is perhaps a warning that the underlying dynamics is more complicated. 
Further study is needed in this area of light-light scalar mesons. 

8.2 Heavy-light quarks 

One of the most promising areas to study scalar mesons on lattice is for heavy-light mesons. 
The scalar meson with cs quantum number is known experimentally [|I|] to be very narrow (it decays 
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Figure 5: Diagrams for «o — > TJ7T and transition strength on a lattice (preliminary results from McNeile et 
al. 



only via the isospin-violating channel D s % or electromagnetically), while the scalar meson with bs 
quantum number is predicted to be similarly narrow from a lattice study of its energy [Eg]. 

The heavy-light scalar meson, which in the limit of a static heavy quark, is expected to be sta- 

: it. d, considered as degenerate). 

►B(0~)7T gives 



ble [26] for bs content and to decay to Bk for bn content (where n ■ 
Evaluating the diagram shown in fig. 6, a lattice estimate of the decay rate of B(0 4 



a width predicted Q27Q as 162(30) MeV. This state has not been observed experimentally yet, but 
the experimental results for the corresponding cn state, D(0 + ), are that the width is 270 ± 50 MeV. 
Although significant 1 /rag effects are expected in the HQET in extrapolating to charm quarks, this 
is indeed a similar magnitude to that predicted for B mesons. It will be interesting so see how the 
lattice prediction of the mass and width of the B(0 + ) fares when experimental results are available. 
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Figure 6: Diagrams for B(0 + ) — > Bn and transition strength on a lattice. 
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9. Do decays matter? 

The previous discussion of decays on a lattice emphasises that qq states do mix with two-body 
states with the same quantum numbers. In the real world, the two-body states are a continuum 
and nearby states have a predominant influence. Since there is a suppression in the amplitude near 
threshold from the factor q L for an L-wave transition, for S-wave transitions (L = 0) the threshold 
will turn on most abruptly and hence will have stronger mixing. 

For bound states there is an influence of nearby many-body states (eg. Ntz on N or 717171 on 
71) which mix to reduce the mass. The nearest such thresholds will be those with pionic channels 
since pions are the lightest mesons. This is the province of low energy effective theories, especially 
Chiral Perturbation Theory which is discussed in other talks. This then provides a reliable guide in 
extrapolating lattice results to the physical light quark masses. 

For unstable states (resonances) the influence of the two-body continuum is less clear since 
the two-body states are both lighter and heavier. In the continuum at large volume, effective field 
theories can again be used to explore this. On a lattice, however, the signal for a particle becomes 
obscured as the quark mass is reduced so that it becomes unstable. Techniques, such as those 
discussed above, are needed to extract the elastic scattering phase shift and hence the mass and 
width. 

For quenched QCD, however, where these two-body states are not coupled (or have the wrong 
sign as in ciq — ► f]7i), then the unstable states will be distorted compared to full QCD. For instance, 
in existing quenched QCD studies, the p will be too heavy since it is not repelled by the heavier %% 
states. Indeed an example of this effect was seen above in the study of p decay including dynamical 



sea quarks, where the p mass decreased [ ]14[ ] when it could couple to %% compared to when it could 
not. 



10. Molecular states? 

Can lattice QCD provide evidence about possible molecular states: hadrons made predomi- 
nantly of two hadrons? 

The prototype is the deuteron: n p bound in a relative S-wave (with some D-wave admixture) 
by 71 exchange. 

There are many states close to two-body thresholds. Since S-wave thresholds are the most 
abrupt, it is usually in this case that the influence of the threshold on the state has been discussed. 
Some of these cases are: 

/o(980) a (980) «-» KK D s (0+) D(0~)K 

B s (0+) «-> B(0~)K X(3872) <-> D*D 
A(1405) «-» KN AT(1535) <-► r]N 

Some of these states (D s (0 + ), B s (0 + )) are stable (in QCD in the isospin conserving limit) 
whereas the rest have other channels open. There is a large literature, stretching over 40 years, 
discussing the consequences of the nearby threshold on these states. One definite implication is 
that isospin breaking is enhanced by mass splittings in thresholds (eg. i^ ^ compared to K + K is 
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8 MeV higher and this induces isospin mixing between the states at 980 MeV). This level of detail 
is not accessible in lattice studies at present, but lattice QCD should be able to address the issue of 
the influence of thresholds on these states. 

The observation of a state near a 2-body threshold implies that there is an attractive interaction 
between the two bodies. But this is a topic like that of whether the chicken or egg was created 
first: an attractive interaction implies and is implied by a nearby state. What can lattice QCD offer 
here? We are in the position of being able to vary the quark masses and this is a very useful tool. 
A two-body threshold will move in general in a different way with changing quark mass than a qq 
state. We can also move the strange and non-strange masses separately and this can be helpful too. 

Another line of investigation is that lattice studies can explore the wavefunction of a state - 
either the Bethe-Salpeter wavefunction or the charge or matter spatial distribution. One can also 
explore the coupling of a state to a 2-body channel, as was discussed above. 

The prototype of a molecular state is the deuteron: it has a tiny binding energy (2.2 MeV) 
and a very extended spatial wave function. Pion exchange between neutron and proton gives a 
mechanism for this long-range attraction. In general it is difficult to reproduce such small binding 
energies in lattice studies. 

Another case where a long-range pion exchange can give binding is in the BB system. Here 



lattice results indicate [ |28| ] the possibility of molecular bound states in some quantum number 
channels which have an attractive interaction from pion exchange, but also the possibility of bound 
multi-quark states which are not described as hadron-hadron but where the two heavy quarks form 
a colour triplet and the light quarks are arranged as in a heavy-light-light baryon. This BB example 
illustrates the rich structure available to multi-quark systems. 



Cases that have been studied on the lattice are the bs and cs scalar mesons |29p. The 
lattice mass values do suggest that these states, treated as Qq states, are bound and are not unduly 
influenced by the BK, DK threshold. A study of the charge distribution of the B(0 + ) gives [30, 31] 
confirmation since the light quark spatial distribution is similar to that of other Qq states. 

For the ciq and fo at 980 MeV, there is limited progress. Quenched studies are inappropriate 



here for the cto since the r\K decay is wrongly treated. In dynamical studies [ ]12| , |32| ], the two body 
channels are of similar energy to the ao (as discussed above), so analysis is unclear. A thorough 
study of the light-light scalar sector is still awaited. 

11. Multi-quark states? 

I told you so: a narrow pentaquark above KN threshold is not possible in QCD. I lived through 
the split Ai and baryonium, which were both narrow features that appeared significant experimen- 
tally but which were evanescent [[33|]. So I have always advised lattice theorists not to rely on the 
experimental evidence for a pentaquark. A pentaquark state around 1540 MeV has a "fall-apart" 
decay mode to KN. This decay does not involve any quark pair production so is expected to be 
unsuppressed - resulting in a very wide decay width (of the order of hundreds of MeV). 

Since the flimsy experimental evidence for the pentaquark is now weakening [jTL [34, 35], it is 



less compelling to study it exhaustively on the lattice. Nevertheless, it is instructive to discuss what 
lattice study could do to make more firm my conviction that "a narrow pentaquark is not possible 
in QCD". The essential issue is the narrowness of the claimed signal, since a broad resonance state 
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would not be unexpected. In a lattice exploration of the pentaquark, an attractive phase shift in 
some KN channel is not sufficient to resolve the issue, the width needs to be evaluated. This needs, 
in principle, 

• 2+1 flavours of sea quark, with light u,d. 

• Operators to create multi-quark states and two-body states 

• Vary spatial size to determine phase shift of two-body interaction. 

Such a comprehensive study is not yet available. Instead most studies have used quenched 
QCD. The KN decay channel is coupled in quenched studies to a pentaquark state, so this suggests 
that a quenched study can be a useful first step, although care must be taken of spurious contribu- 
tions arising from states such as KNf] which will be relatively light and have unphysical behaviour. 
For realistic quark masses, there will be two body states on the lattice: their energy shifts then, 
following Luscher, give the phase shift. By exploring this phase shift versus energy, the width 
of any resonance can then be extracted. In practice, for a very narrow resonance, the spectrum 
will look more like a one-body state (the resonance) and a collection of two-body states except for 
mixing near the avoided-level crossings. As the lattice volume is varied, the two body levels with 
non-zero momentum will move (since momentum is 2%njL ). Thus the volume dependence is a 
useful diagnostic. It is difficult, however, to determine accurately the many energy levels expected 
in a relatively large volume. An additional indicator is that the weight of the one-body and two- 
body contributions can be volume dependent: basically because the contribution from a two body 
state to a local-local correlator will be dominated by the lowest relative momentum which has a 
contribution which behaves as 1/L compared to the contribution from a one-body state. 

Using these criteria, lattice groups have explored the S- and P-wave KN system. A summary 
of lattice results at LAT04 [ |36| ] showed most groups reproducing a pentaquark state. Many more 
results have been presented since and at this conference. Here I comment on some of the more 



comprehensive studies. The Kentucky group [37] use relatively light quarks and conclude no evi 



dence for any pentaquark. One group which claimed lattice evidence for a pentaquark signal [38] 
has subsequently withdrawn their claim [J3"9[]. The analysis of quenched lattices at quark masses 
heavier than physical is still somewhat subjective and lattice evidence for a pentaquark is still being 



claimed by some groups, e.g. [40]. 



12. Conclusions 



Hadronic physics involves unstable states. 

Don't put your head in the sand: study these with lattice techniques. 

We need to study interactions as well as spectra: solid gold - not just gold-plated. 
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